Green's Formula 



In this note we state Green's formula and look at some examples. We 
will prove it in the next note. 

1. Green's Formula 

Suppose that we have a linear time invariant system with rest IC. 

P(D)y = f(t), y(f) = for f < (1) 

• As in previous sessions, we will consider /(f) to be the input to this 
system. Everything we say will also hold for systems like 

T + kT = kT e with input T e and mi + bx + kx = by with input y. 

• In this context, where we don't consider functions for f < 0, the initial 
conditions mean that j/(f) and all its derivatives are at f = 0~ . 

• P(D) is a polynomial differential operator. Although it can be of any 
order, recall that we developed the second order case extensively in 
the last unit, where it was often written as 

P(D)y = my + by + ky. 

Suppose further that w(t) is the unit impulse response for (1). That 
is, w(t) satisfies P(D)w = S(t), with rest IC. Then, for any input /(f) the 
solution to equation (1) is given by Green's formula 

y(t) = (J*w)(t) = f f{x)w{t-x)dx. (2) 

This is a wonderful formula! It tells us the response to any input once we 
know the unit impulse response. Furthermore, it gives us that response as 
an integral which can be computed numerically if necessary. For many 
physical systems the impulse response can be measured directly or de- 
duced from measurements. So, Green's formula gives us a method for 
predicting the system's response to any input. 

2. Unit Impulse Response = Weight Function 

The unit impulse response is also called the weight function. We will 
use the terms interchangeably. If we think of an integral as a 'sum' then 
Green's formula shows the solution y(t) to (1) is given as a weighted sum 
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of the small bits of input, /(t) dr from before time t. Each piece is weighted 
by w(t — t). 

Before proceeding, let us recall the definition of the unit impulse re- 
sponse. The weight w(t) is the unique solution to the IVP 

P(D)y = 5(t) with rest IC (3) 

In the previous session we learned how to rewrite (3) as a homegeneous 
equation. We will only restate this for second order equations. The weight 
function for the system 

mx + bx + kx = /(f) 
is for t < and the solution to 

mx + b'x + kx = 0, x(0) = 0, x(0) = 1/m 

for t > 0. 

3. Examples 

We now out Green's formula (2) in a couple of cases where we can check 
it by finding the particular solution y v by another method. 

Example 1. Find the particular solution given by (2) to 

y + y = A, y(0) = 0, y (0) = 0, where A is a constant. 

Solution. The unit impulse response is w(t) = sinf. Therefore for t > 0, 
we have 



Vv 



(t)= [ Asm(t - t) dT = Acos(£ - r) 
Jo 



t 

= A{1 -cost). 

J o 



We check this by another method: The exponential response formula or 
the method of undetermined coefficients produces the particular solution 
y v = A. Adding in the homogeneous solution we get the general solution 
to the DE is 

y = A + c\ cos t + C2 sin t. 

You can easily compute that the rest initial conditions are matched by 
y = A — A cos t, as found by Green's formula. 
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Example 2. Find the particular solution for t > given by (2) to 

y"+y = f{t) = 



1 for < t < 7i 
elsewhere 



Solution. Here the method of Example 1 leads to two cases: < t < n and 
t > n: 



Vv= I /(t) sin(t - t) rfr = 

J 



sm{t — t) dr = cos(t — t) 
Jq sin(f — t) dr = cos(t — t) 

We leave it to ther reader to check this by our earlier methods 



1 — cos t, for < t < n; 



= —2 cos t, for t > n . 



J 
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